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Abstract. We prove three conjectures concerning the evaluation of de- 
terminants, which are related to the counting of plane partitions and 
rhombus tilings. One of them has been posed by George Andrews in 
1980, the other two are by Guoce Xin and Christian Krattenthaler. Our 
proofs employ computer algebra methods, namely the holonomic ansatz 
proposed by Doron Zeilberger and variations thereof. These variations 
make Zeilberger's original approach even more powerful and allow for 
addressing a wider variety of determinants. Finally we present, as a chal- 
lenge problem, a conjecture about a closed form evaluation of Andrews's 
determinant. 
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1. Introduction 

The concept of determinants evolved as early as 1545 when Girolamo Car- 
dano tried to solve systems of linear equations. The mathematics community 
slowly realized the importance of determinants; we had to wait for more than 
200 years before someone formally defined the term "determinant" . It was 
first introduced by Carl Friedrich Gaufi in his Disquisitiones Arithmeticae in 
1801. Determinants possess many nice properties and formulas such as multi- 
plicativity, invariance under row operations, Cramer's rule, etc. Every student 
nowadays learns how to compute the determinant of a specific given matrix, 
say with fixed dimension and containing numeric quantities as entries. On 
the other hand, there are lots of matrices with symbolic entries that have a 
nice closed-form formula for arbitrary dimension. The first example in the 
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history of mathematics and still the most prominent one is the Vandermonde 
matrix. 

Starting in the mid 1970's, the importance of evaluating determinants 
became even more significant when people related counting problems from 
combinatorics to the evaluation of certain determinants. Evaluating determi- 
nants of matrices whose dimension varies according to a parameter is a deli- 
cate problem but at the same time a very important one. Around the same 
time the field of computer algebra emerged. Doron Zeilberger was amongst 
the first mathematicians to realize the importance of computer algebra algo- 
rithms for combinatorial problems, special function identities, symbolic sum- 
mation, and many more. In his paper |17j he built the bridge between these 
two topics, namely symbolic determinant evaluation and computer algebra, 
and since then his holonomic ansatz has been successfully applied to many 
problems related to the evaluation of determinants. The most prominent one 
is probably a notorious conjecture from enumerative combinatorics, the so- 
called g-TSPP conjecture, which was the only remaining open problem from 
the famous list [H] by Richard Stanley (it also appeared in [T^]) until it was 
recently proved |10| using Zeilberger's holonomic ansatz. 

In this paper, we solve some of the problems that are listed in Christian 
Krattenthaler's complement [TH] to his celebrated essay [H] (the attentive 
reader may already have observed that our title is an allusion to this refer- 
ence). At the same time we show that the holonomic ansatz can be modified 
in various ways in order to apply it to particular problems that could not be 
addressed with the original method. 



2. Zeilberger's Holonomic Ansatz 

For sake of self-containedness we recall briefly the original holonomic ansatz 
for determinant evaluations as it was proposed by Zeilberger |17l . Its steps are 
completely automatic and produce a rigorous proof — provided that they can 
be successfully carried out in the example at hand. In particular, the approach 
relies on the existence of a "nice" description for an auxiliary function (it 
appears as c„ j below); if such a description does not exist then the holonomic 
ansatz fails. That's why we call it an "approach" or an "ansatz" , rather than 
an algorithm. 

Generally speaking, Zeilberger's holonomic ansatz addresses determi- 
nant evaluations of the type 

AeiAn — det (a^j ) = &„ (n > 1) 

where the entries of the n x n matrix An and the (conjectured) evalu- 
ation bn (where 6„ 7^ is required for all n ^ 1) are explicitly given. The 
underlying principle is an induction argument on n. The base case ai.i = 61 is 
easily checked. Now assume that the determinant evaluation has been proven 
for n — 1. In particular, it follows that det74„_i is nonzero by the general 
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assumption on fe„. Hence the rows of An-i are linearly independent and thus 
the linear system 



/ ai,i 



On- 1.1 




ai 



^n— l.n— 1 ^n — l,\ 
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has a unique solution. In the well-known Laplace expansion formula (here 
with respect to the last row) 



det (a,j) 



(-l)"+^M„^,a„,, 



the expression (— l)"+^M„j- is called the (n, j)-cofactor of An- The minor 
Mnj is the determinant of the matrix obtained by removing the n-th row 
and the j-th column. The above linear system has been constructed in such 
a way that the entry c„j- in its solution is precisely the (n, j)-cofactor of An 
divided by its (n, n)-cofactor (which is just the (n — l)-determinant). This fact 
can easily be seen by Cramer's rule, i.e., by considering the matrix A^l^ that 
we obtain from An by replacing the last row by the z-th row. For 1 ^ i < n 
the Laplace expansion of An^ corresponds exactly to the i-th equation in our 
linear system. Now the determinant of An is given by 

n 

^n — 1 ^ ^ ^n. j^n.j ■ 
J = l 

To complete the induction step it remains to show that this expression is 
equal to bn- 

The problem is that we cannot expect to obtain a closed-form expres- 
sion for the quantity Cn.j (otherwise we certainly would be able to derive a 
closed form for the determinant and we were done). Instead, we will guess 
an implicit, recursive definition for a bivariate sequence c„j and then prove 
that it satisfies 



n 

Cn,jO'i,j = 



J = l 



(n ^ 1), 
(1 s$ i < n), 

{n ^ 1). 



(1) 
(2) 

(3) 



From the first two identities which correspond to the linear system given 
above, it follows that our guessed c„j is indeed the normalized (n, j)-cofactor. 
Identity ^ then certifies that the determinant evaluates to 5„. Hence the 
sequence Cnj plays the role of a certificate for the determinant evaluation. 

Now what kind of implicit definition for c„j could we think of? Of 
course, there was a good reason for Zcilberger to name his approach the 
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"holonomic ansatz". That is to say because he had the class of holonomic 
functions (or better: sequences) in mind when he formulated his approach. 
In short, this class consists of multi-dimensional sequences that satisfy linear 
recurrence relations with polynomial coefficients, such that the sequence is 
uniquely determined by specifying finitely many initial values (we omit some 
additional technical conditions here). What makes the use of this class of se- 
quences convenient is the fact that it is not only closed under the basic arith- 
metic operations (addition, multiplication), but also under more advanced 
operations such as specialization, diagonalization, and definite summation; 
all these are needed to prove ([T]), ([2]), and ([3]). Moreover, there exist known 
algorithms for performing zero tests and the previously mentioned opera- 
tions. For more details on holonomic functions and related algorithms see, 
e.g., [T3 m H]). Thus if the matrix entries Uij are holonomic and if luckily 
the auxiliary function c„ j turns out to be holonomic, then the approach will 
succeed to produce a recurrence for the quotient 6„/&„_i. 



3. An Old Problem by George Andrews 

In the context of enumerating certain classes of plane partitions, namely cycli- 
cally symmetric ones and descending ones, George Andrews [5] encountered 
an intriguing determinant which he posed as a challenging problem. In Chris- 
tian Krattenthaler's survey [13j . it appears as Problem 34- This conjecture 
does not even deal with a closed-form evaluation, but it is "only" about the 
quotient of two consecutive determinants; a situation that strongly suggests 
to employ Zeilberger's holonomic ansatz! 

Theorem 1. Let the determinant Di{n) be defined by 

''fi + i+j -2 



Di{n) :— det 5i , 

l^i,j<n \ V ] 

where fi is an indeterminate and Sij is the Kronecker delta function. Then 
the following relation holds: 



= (-1) 



Proof. By looking at the first few evaluations of Di (n) (for 1 ^ n ^ 8 they are 
explicitly displayed in [H]), we see that only the quotient Di{2n) / Di{2n — 1) 
is nice, but not Di{2n+1) / Di{2n) . The reason is the occurrence of irreducible 
nonlinear factors that change every two steps, i.e., Di{2n) and Di{2n — 1) 
share the same "ugly" factor (thus their quotient is nice), but in Di (2r7,4-l) the 
"ugly" part will be different (and therefore the quotient Di{2n + 1)/ Di{2n) 
does not factor nicely). 

We first tried Zeilberger's original approach on the determinant Di{n). 
But we didn't even succeed to guess the recurrences for c„ j- as they either are 
extraordinarily large or do not exist at all. Moreover, we have good reasons to 



Advanced Computer Algebra for Determinants 



5 



believe that the quotient Di{2n + 1) / Di{2n) is not holonomic (see Section 15]), 
in which case we know a priori that the approach cannot succeed. 

Therefore we have to come up with a variation of Zeilberger's approach, 
that pays attention to even n only. This means that we consider the normal- 
ized cofactors c„j only for matrices of even size. The new identities, to be 
verified, will be 

{n ^ 1), (la) 
(1 i < 2n), (2a) 

(n ^ 1). (3a) 

In order to come up with an appropriate guess for the yet unknown 
function c„j-, we compute the normalized cofactors for all even-size matrices 
up to dimension 30. This gives a 15 x 30 array with values in (Q(/x) that is 
used for guessing linear recurrences for c„j . For this step, Manuel Kauers's 
Mathematica package Guess has been employed, see [7| for more details. To 
give the reader an impression of what the output looks like, we display the 
results here in truncated form: whenever the abbreviation {k terms) appears, 
it indicates that this polynomial cannot be factored into smaller pieces, and 
that for better readability it is not displayed in full size here. 



Cn,2n — 1 



'^njO-iJ — 

2n 



E 



2n{j + l)i2n ~ l){2j + - 2?i)(j - 2n + 1) 

X (/i + An - 5)(/i + 4n - 3)(j + ^ + 2n - l)c„j = 
j{j +n- l){2j +^Ji- l)(j -2n + 3)(m + 4n - 3) (4) 

X {f + 2jV + ■ • • (24 terms) + 12)c„_i,j+i - 
U + l)(j + M + 2n - 3)(2jV + 8j®n +...{92 terms) - 210yun)c„_ij 

(j - l)(j + M - 3)(j + 2){2j + 2n){j + fi + 2n- l)c„,, = 

j{j + fi- 3)(4j4 + 8ffi + ... (26 terms) + 16)c„j_i - (5) 
jU - 1)(J + - 2)(2j 2){j -2n- 2){j + ^i + 2n- 3)c„,j_2 

When translated into operators, the two recurrences and (O consti- 
tute a Grobner basis of the left ideal that they generate in the corresponding 
operator algebra. As a consequence it follows that they are compatible; this 
means that starting from some given initial values, these recurrences will 
always produce the same value for a particular Cngjo, independently from 
the order in which they are applied. The support of the recurrences suggests 
that fixing the initial values ci.i ~ — /i/2 and ci_2 = 1 is sufficient: ([5]) would 
produce cij for all j > 2, and then H]) could be used to obtain the full array 
of values (c„j)„j^i. 
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Unfortunately it is not that easy, since there are two disturbing pheno- 
mena. The first is the factor j — 2n that appears in both leading coefficients. 
Hence for computing c„_2n none of the recurrences is applicable and we are 
stuck. Note that the factor j — 2n+ 1 in the leading coefficient of ([4]) is not 
a problem since we can use ([5]) instead. We overcome the problem with c„_2n 
by finding a recurrence of the form 

P2C„+2.,j+4+PlC„+l.j+2+PQC„.j ^0, P0,?3l,P2 e (Q[n,j,/X] (6) 

in the ideal generated by ^ and (O. The coefficient p2 does not vanish for 
j = 2n, and thus, together with the additional initial value C2,4 = 1, the 
recurrence ([6|) allows to compute the values Cn,2n- 

The second unpleasant phenomenon is the free parameter fi which can 
cause the same effect for certain choices. For example, if /i = —6 then we can- 
not compute C2,3 from the given initial values since both leading coefficients 
vanish simultaneously by virtue of the factors j — 2n -\- 1 and j -\- fj, + 2n — 1. 
We handle it by restricting the parameter /i to the real numbers > /ig for a 
certain /ig S K, and by showing that all our calculations are sound under this 
assumption (in most cases we use /i > or > 2). But since the determinant 
for every n G N is a polynomial in fj, we can extend our result afterwards 
to all /i e (D. Alternatively, one can argue that is a formal parameter and 
therefore expressions like fj,+6 are not zero as they are considered as elements 
in the polynomial ring (D[/i]. 

We have now prepared the stage for executing the main part of Zeil- 
berger's approach. Identity ([Ta|) is easily shown with help of recurrence 
Substituting j 2n produces a recurrence for the entries c„^2n, which to- 
gether with the initial values implies that Cn,2n = 1 for all n. The identi- 
ties (Pa)) and ([5a|) are proven automatically as well, since they are standard 
applications of holonomic closure properties and summation techniques (cre- 
ative telescoping). For these tasks we have used the first author's Mathe- 
matica package HolonomicFunctions 9 . The interested reader is invited to 
study our computations in detail, by downloading the electronic supplemen- 
tary material from our webpage [IT]. □ 

4. Interlude: The Double- Step Variant 

In this section we propose a variant of Zeilberger's holonomic ansatz, de- 
scribed in Section O that enlarges the class of determinants which can be 
treated by this kind of ansatz. The condition 6„ ^ for all n ^ I im- 
poses already some restriction. For example, when studying a Pfaffian Pf {A) 
for some skew-symmetric matrix A, one could be tempted to apply Zeil- 
berger's approach to the determinant of A; recall that Pi{A)^ = det{A). The 
problem then is that det{A) is zero whenever the dimension of A is odd. 
Hence one would like to study the quotient 6„/&„-2 instead of the forbidden 
bn/bn~i'i as in the previous section, n has to be restricted to the even integers. 
This dilemma can be solved by modifying Zeilberger's ansatz subject to the 
Laplace expansion of Pfaffians, see [B]. 
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On the other hand there are determinants which do factor nicely for 
even dimensions but not for odd ones. Also here, we expect the quotient 
bn/bn-2 to be nice, whereas the expression 5„/&„_i might not even satisfy a 
linear recurrence and hence could not be handled by Zeilberger's holonomic 
ansatz at all. Similarly when the closed form 6„ is different for even and 
odd n: while here the original approach could probably work in principle, 
one may not succeed because of the computational complexity that is caused 
by studying the quotient &„/&„_i, which is expected to be more complicated 
than bn/bn-2- See Theorems [2] and [5] below for such examples, which actually 
have motivated us to propose the following variant. 

As we announced we now generalize Zeilberger's method in order to 
produce a recurrence for the quotient of determinants whose dimensions differ 
by two. As before, let M = (aij)i^ij-^„ and let bn denote the conjectured 
evaluation of det(A/), which for all n in question has to be nonzero. We pull 
out of the hat two discrete functions c'^ ^ and c" j and verify the following 
identities: 



0, (lb) 



n 



E = E =0 (2b) 

. ^^^^ 

Then the determinant evaluation follows as a consequence, using a similar 
induction argument as in Section [51 

Let's try to give the motivation for this approach which also explains 
why it works. The idea is based on the formula for the determinant of a block 
matrix: 

det(M) = det f ^{^ f,^ ) = det(Afi) det(M4 - MgAff^Ma). 

We want to divide the matrix M into blocks such that M4 is a 2 x 2 matrix. 
Let C — (C, C") denote the (n — 2) x 2 matrix whose first column is C" = 
and whose second column is C" — (c"j)i^j^n-2- With this 
notation, the two equations (pb)) can be written as 

' C" C" 

(Ml, M2) I 1 1 = (Afi, M2) f f I = MiC + M2 = 0, 



1 



where the conditions of Equation (llb[) have been employed to constitute the 
identity matrix 12- Now by the induction hypothesis we may assume that 
deti^i,j^n-2i(ii,j) = det(Afi) equals 6„_2, which by our general assumption 
is nonzero. Thus the above system determines C uniquely and we can write 
C = -M^^Nh. 
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Finally, we obtain the missing part from the block matrix formula that 
gives us the quotient det(Af )/ det(A/i), as the determinant of the 2x2 matrix 

(M3, M4) ) = + M4 ^ Ah - KhM^^Ah. 

This is exactly what is expressed in Equation (j3bp . with all matrix multipli- 
cations explicitly written out. 

Even though it turned out during our research that there is no need to 
apply the double-step method in the present context, we decided to keep it, 
as we are sure that it will be useful for future use. Let us also remark that 
it is straight-forward to generalize this idea in order to produce a triple-step 
method, etc. But since the identities to be verified become more and more 
complicated, we don't believe that these large-step methods are relevant in 
practice. 



5. Small Change, Big Impact 

The next determinants we want to study appear as Conjecture 35 and Con- 
jecture 36 in |13) . Christian Krattenthaler (private communication) describes 
the story of how they were raised: 

I wrote this article during a stay at the Mittag-Leffler Institut. 
Guoce Xin was also there, as well as Alain Lascoux. They followed 
the progress on the article with interest. So, it was Guoce Xin, 
who had been looking at similar determinants at the time (in the 
course of his work with Ira Gessel on his big article on determi- 
nants and path counting [5 ), who told me what became Conjecture 
35. I made some experiments and discovered Conjecture 36. Alain 
Lascoux saw all this, and he came up with Conjecture 37. 

Guoce Xin's observation was that a certain matrix, very similar to the 
one of Section [31 has a determinant that factors completely; the only change 
is the sign of the Kronecker delta 5i^j. But still, the evaluation is given as a 
case distinction between even and odd dimensions of the matrix. 



Theorem 2. Let fi be an indeterminate and n be a nonnegative integer. Then 
the determinant 
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if n is even, and it equals 

^ ' V2^^ ^')in+i)/2[ l\ (2z)!(2z + 2)!^ 



L(«+1)/4J 

n i(/^ + 6i-l) (i(-/i-3n + 6i-3) 

V^' J (n-4i+l)/2 \^ /(ri-4i+3)/2 



if n is odd. 

Proof. We could try to solve this problem directly, either by Zeilberger's 
original ansatz or in the way we did in Section [31 However, this does not 
work in practice as the computations become too large (in the second case we 
were even able to guess the recurrences for c„j , but their size destroyed any 
hope to prove ([2a|) and ([3a|) ). Next, we could give the double-step method 
a try, which we described in Section 21 We succeeded to make the proof 
for some concrete integer fi, but again it seems intractable for symbolic /i. 
Instead we make a small detour and break this determinant into pieces in 
order to make the calculations smaller. To obtain the desired result, we put 
these pieces together by the Desnanot-Jacobi adjoint matrix theorem, also 
known as Dodgson's rule; this formula gave rise to a celebrated algorithmic 
determinant evaluation method as well [16, 1 , but that approach is different 
from our usage of Dodgson's rule. Let us introduce the following notation 

^I + i+j -2^ 

j 



bn{I, J) := bnil, J, := ^i'^^^i+i [ ^'^^J 

SO that our determinant ^ is denoted by 6„(1,1). In this notation the 
Desnanot-Jacobi identity is stated as follows: 

&„(0,0)6„_2(1,1) = 6„-i(0,0)fo„_i(l,l)-&„_i(0,l)fe„_i(l,0). (8) 

After substituting n — > 2n + 2 and n ^ 2n + 1 in ^ and using the fact 
^2n(0, 0) = — fo2n-i(l, 1) (to be showu later) gives the equations 

Wi(0,l)62«+i(l,0) 



52„+l(l, 1) 
fo2„-l(l,l) 



fo2„(l,l)+ Wi(o,o) 

-&2„(0,1)62„(1,0) 



fo2„(l,l) + 62„+l(0,0) 

from which the desired quotient can be obtained: 

62„+i(l,l) Wi(0,l) 62„+i(l,0) 



(9) 



^'2„-l(l,l) 62„(0,1) 62«(1,0) 

Similarly, we substitute n ^ 2n + 1 and n 2n into ([8]) and use the fact 
62n-i(0,0) — (again, to be shown later), to derive the quotient of even 
determinants 

^>2n(l,l) ^ &2»(0,1) 62n(l,0) . . 

62„-2(l,l) 62„-l(0,l) ■ 62n-l(l,0)' ^ ' 

We now use the first variation of Zeilberger's ansatz (see Section [3]) to derive 
recurrences for the quotients 62n-i-i(0, l)/62n(0, 1), etc. which appear on the 
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right hand sides of © and pUj) . Since the arguments closely follow the lines 
of the proof of Theorem [T] we don't detail further this part and refer to the 
electronic material jTT]. Although for our purposes it would suffice to work 
with these recurrences, we succeed in solving them in closed form: 

^ , , _ 62„+i(0, 1) _ 2 (f + 2n)„^,(/i + 2n- 1^^ 

62„+i(l,0) 2(f + 2n)„^^(/i + 2n + l)„_^ 
b2n(l,0) 2(f + 2n-l)„_^(/z + 2n+l)„_^ 

- ra^ ^ — — ■ '''^ 

These quotients immediately give closed form evaluations of the correspond- 
ing determinants (see also Theorems [3] and E]) . It remains to justify the as- 
sumptions 62,1(0,0) = —62,1-1(1,1) and 62,1-1(0,0) = that were used to 
derive © and pO)) . 

In order to evaluate the quotient 62,1(0, 0)/62„-i(l, 1) we need to modify 
the method presented in Section |31 we apply Laplace expansion with respect 
to the first row of the matrix, instead of the n-th row, and we normalize 
the auxiliary function c„j in such a way that Cnfi = 1. If we come up with a 
recursive description of some function c„ j and are able to verify the following 
identities, then we are done: 



c„,o = 1 {n> 1), (Ic) 

2n-l 

^ Cn.ja^.j =0 (0 < i 2n - 1), (2c) 

62n(0,0) 1 . r.^ 

2^ c„,jao,j = 7 77-rT = -1 {n> 1). (3c) 

~^ 62„-l(,l,lj 

As before, the details can be found in [TT]. 

Last but not least we have to argue that 62n-i(0, 0) vanishes. Christian 
Krattenthaler kindly pointed us to Theorem 11 in ^ (see also Theorem 35 
in |12j ) which contains this statement. Anyway we have produced an alterna- 
tive, computerized proof: actually it is very simple, since we just have to come 
up with a guess for the coefficients of a linear combination of the columns 
(or rows) that gives 0, and then prove that our guess does the job. Hence we 
find a recursive description of a function Cn,j for n ^ 1 and ^ j ^ 2n — 2, 
such that 

Cn,0 • + • ■ • -|- C„ 2n-2 



\a2n- 



2,oy 



^ ao,2,i-2 ^ 




3 


\a2,i-2,2n~2/ 
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and such that there is an index j for which Cnj ^ 0. For guessing, we compute 
Cn.j for all n up to some bound and normalize them. Luckily the nullspace 
of the above system has always dimension 1, otherwise it would be trickier 
to find a suitable linear combination (however, in our problem this is no 
surprise since we are finally aiming at proving that the minor Mq^q of this 
matrix evaluates to some nonzero expression) . So we are successful in guessing 
the recurrences of c„j and use them to prove 

C„,2n~2 = 1 {n^ I) 

2n-2 

Cnja.j =0 (0 ^ i ^ 2n - 2). 

This concludes the proof of Theorem [5] □ 

Since the evaluations of the determinants 6„(0, 1) and &„(1,0) are in- 
teresting results in their own right, but are somehow hidden in the proof of 
Theorem [21 we are going to state them explicitly here. 

Theorem 3. Let fi be an indeterminate and n be a nonnegative integer. Let 
Qi{n) and Q2{n) be defined as in (|11|) and ()12|) . respectively. Then the de- 
terminant 

6„(0,1)= det (^S,_,,^ + (^^'^J'^ 



/ n/2-1 \ / n/2 ^ 

n n^^w 

\ fc=0 / \fc=l / 

if n is even, and it equals 

(n-l)/2 

n (Qi(fc)Q2(fc) 

fe=l 

if n is odd. 

Proof. It is analogous to the proof of Theorem[T]and can be found in [TI] . □ 

Theorem 4. Let fi be an indeterminate and n be a nonnegative integer. Let 
Qsin) and Qi{n) be defined as in (jl3p and (|14p . respectively. Then the de- 
terminant 

6„(1,0)= det (^S.^^_,+ (^' + '+^-^ 



n/2-1 \ / n/2 ^ 

n Q3(fc) n 

k=0 ) y 
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if n is even, and it equals 

(n-l)/2 

n {Qsimm 

k=l 

if n is odd. 

Proof. It is analogous to the proof of Theorem [T] and can be found in [H] . □ 

By replacing j by j + 1 at the bottom of the binomial coefficient in the 
entries of ([7]), we arrive at our last determinant; it has been discovered by 
Christian Krattenthaler and appears as Conjecture 36 in his paper [13 (note 
that the formula there is erroneous, as one of the product quantors is missing 
so that the corresponding factor slided into the previous product). Also this 
problem has its own combinatorial interpretation in terms of counting certain 
rhombus tilings. 

Theorem 5. Let fi be an indeterminate. For any odd nonnegative integer n 
there holds 

det u.+f'^'-'tr'' 

KiJ^n \ V J + 1 

C_n("-l)/2 2(»-l)("+5)/4 (.. I 1^ ^^^^ ^ 2))(«+i)/2 

^ ^ ^ (^(- + 1))! 

L(«+3)/4J _ 



(ri-4i+l)/2 





n {li-i^-^ 
1=1 




Proof. Because of the similarity of this determinant with ([7]) , we are able to 
relate these two problems via shifting the starting points: 



det I —6ij 



l^i,i<2n-l \ \ i + 1 

'{n-2)+i+j -2 



det —Sii 

2^i,j^2n V V j 

By using the notation from above, the determinant of Theorem [5] is denoted 
by 62n-i(2, 2,11 — 2). Analogously to P^ - (l5c)) . we apply a variation of Zeil- 
berger's approach to derive a recurrence for 

_ &2n(l, 

'^"^''^"62„-i(2,2,/z)- 

The result is qn+i{fJ-) — qnip-) — which reveals that the quotient g„(/i) is 
constant. Together with the initial value = — 4/(/i + 3) and the fact 

that 62n(l, 1, m) is already known from Theorem[21 we get the desired result. 
Once again, we refer to IJA, for the details of the computations. □ 
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As mentioned above, Alain Lascoux found that the more general deter- 
minant 

det 

factors completely for odd natural numbers n and r, and its complicated 
evaluation, which was figured out by Christian Krattenthaler, appears as 
Conjecture 37 in [13] . We remark that the formula given there holds only for 
r n; otherwise the Kronecker delta does not show up in the matrix and the 
evaluation is much simpler. We cannot attack this determinant directly with 
Zeilberger's ansatz since the matrix entries do not evaluate to polynomials 
in /i for concrete integers i and j, as long as r is kept symbolically. Therefore 
the guessing for c„j will not work. A different strategy would consist in 
finding some connection between the cases r and r + 2; then induction on r 
would provide a proof, using Theorem[2]as the base case r = 1. Unfortunately 
we were not able to achieve this goal. 



6. A Challenge Problem 

We want to conclude our article with a challenge problem for the next gener- 
ation of computer algebra tools. In Section[3]we have only proven a statement 
about the quotient of two consecutive determinants (Theorem [T]). But so far 
nobody has come up with a closed form for the determinant Di (n) . We now 
present a conjectured closed form, which, however, we are unable to prove 
with the methods described in the present paper. We have already remarked 
in Section [3] that the quotient Di{n) / Di{n — 1) most probably is not holo- 
nomic; in that case Zeilberger's holonomic ansatz and our variations of it are 
not applicable. 



Conjecture 6. Let fi he an indeterminate and let the sequences C(n), F{n), 
and G{n) he defined as follows: 

1=1 

/LiU("-l)J-2j 2\H^+2)\\ 

Ein)^{p + l)^i n (m + 2z + 6) M 



14 



Christoph Koutschan and Thotsaporn "Aek" Thanatipanonda 



/ Li("-i)J \ 

Fm{n) =i II {fx + 2i + n + mf-^'-^ I 

X JJ (/i - 2i + 2n - 2m + l)i-2»-'" 



F{n) 



Pi{n) 
P2{n) 
G{n) 



E{n)Fo{n) 



if n is even 



i(»-5) 



E{n)Fi (n) JJ (/x + 2i + 2n - 1) if n is odd 



i=l 



T{k) = 55296fc6 + 41472(^ - l)k!^ + 384(30^*^ - QQii + 53)fc^ 
+ 96(At - l)(15/x2 - 42^ + 61)P 
+ 4(19/i^ - 122;u^ + 419;u^ - 544/i + 72)A;^ 
+ (m - 1)(m^ - 14/i^ + lOl/i^ - 160m - 84)A; 
+ 2(m-3)(m-2)(/x-1)(m+1) 

= ^ (26'=(/. + 8fc - I) + 5)),,_3 (i(M + 4/c + 2)),_2 

X + 4fc + 2)),_,,_2 T(fc))/((2fc)! + 6fc - 3))3,^, ) 
S,{n) = ^ (26'=(M + Sfc + 3) (i)^, + 5)),,_, + 4fc + 4)),_, 

+ 4k + 4)),„_,,_, r(A: + i))/ ((2fc + 1)! + 6fc + 1)) 



k=l 



3k+5 



= 2 



3n- 



_^(i(^ + 6n-3))3„_2 /(I (^ + 2)) 



(i(M + 5)),„_3 



(/^ + 3) 



2n-2 _^ Mm - l)-S'l(n) 



^ + 6n + 1))3„_^ / (m + 14) + 4)) 



213 J 

2n-2 , - 1)-S'2(n) 



(5(M + 5)),„_, V 

Pi (n + 1)) if n is odd 
P2(§) if n is even 

Then for every positive integer n we have 
^(i + i + j-2^ 



(/. + 7)(m + 9) 



29 



det I Si 



J 



C(n)P(n)G([i(n+l)j; 



Let us add a few remarks on our conjectured closed form. The elements 
of the sequence C(n) are rational numbers of the form 1/k where k is an 
integer. The sequences F{n) and G{n) consist of monic polynomials in 
with integer coefficients. The F{n) factor completely into linear factors of 
the form {/j. + k) where € N, and thus have positive coefficients. The G{n) 
have positive coefficients as well, but turn out to be mostly irreducible; the 
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only counterexample that we found is G(4) — {fi + 34) (/i'^ + + 954/i + 
5928). They correspond to the "ugly factors" mentioned in Section[31 For the 
convenience of the reader, we provide the Mathematica code for all quantities 
introduced in Conjecture [6] in the supplementary material [VL^- 

In order to come up with this complicated conjecture, we computed the 
determinants Di (n) for 1 ^ n ^ 295 which gave us the first 148 polynomials 
of the sequence G{n). These data enabled us to guess recurrences for the 
subsequences Pi{n) and P2{n); the recurrences by the way are used in [TT] to 
provide a fast procedure for computing Di(n). The Maple package LREtools 
was able to find "closed form" solutions which, after lots of automatic and 
manual simplifications, became the formulae for T{k), Si{n), and Pi{n). 
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